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Q ABSTRACT 

We present a re-analysis of archival HST/NICMOS transmission spectroscopy 

! , of three exoplanet systems; HD 189733, GJ-436 and XO-1. Detections of several 

p I molecules, including H2O, CH4 and CO2, have been claimed for HD 189733 and XO-1, 

but similarly sized features are attributed to systematic noise for GJ-436. The data 
consist of time-series grism spectra covering a planetary transit. After extracting light 
curves in independent wavelength channels, we use a linear decorrelation technique 
account for instrumental systematics (which is becoming standard in the field), and 
measure the planet-to-star radius ratio as a function of wavelength. We use a resid- 
ual permutation algorithm to calculate the uncertainties, in an effort to evaluate the 
effects of systematic noise on the resulting transmission spectra. For HD 189733, the 
uncertainties in the transmission spectrum are significantly larger than those previ- 
ously reported. We also find the transmission spectrum is considerably altered when 
y—( using different out-of-transit orbits to remove the systematics, when some paramc- 

^ ters are left out of the decorrelation procedure, or when we perform the decorrelation 

C*") with quadratic functions rather than linear functions. Given that there is no physical 

reason to believe the baseline flux should be modelled as a linear function of any par- 
C" ticular set of parameters, we interpret this as evidence that the linear decorrelation 

technique is not a robust method to remove systematic effects from the light curves 
f*^ for each wavelength channel. For XO-1, the parameters measured to decorrelate the 

y—i light curves would require extrapolation to the in-transit orbit to remove the system- 

atics, and we cannot reproduce the previously reported results. We conclude that the 
resulting NICMOS transmission spectra are too dependent on the method used to re- 
move systematics to be considered robust detections of molecular species in planetary 
•rH atmospheres, although the presence of these molecules is not ruled out. 

5_i Key words: methods: data analysis, stars: individual (HD 189733), stars: individual 

(GJ-436), stars: individual (XO-1), planetary systems, techniques: spectroscopic 
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1 INTRODUCTION 

The number of known exoplanets is increasing rapidly, re- 
vealing a diverse range of systems with vastly different prop- 
erties. Transiting exoplanets offer a unique opportunity to 
study the structure and composition of planets other than 
those in our own solar system, as transit light curves allow 
a measurement of the radius of the planet, and the comple- 
mentary radial velocity technique measures the mass. From 

* E-mail: Neale.Gibson@astro.ox.ac.uk 



the derived density, the bulk composition of the planet may 
be inferred. 

Transiting planets also provide the opportunity to mea- 
sure the composition of planets' atmospheres. Typically, 
transit light curves are modelled assuming the planet as an 
opaque disk, whose size is defined by the altitude at which 
the atmosphere becomes opaque to starlight. However, the 
optical depth in the atmosphere is wavelength dependent, 
being sensitive to atomic and molecular absorption. There- 
fore, the transit depth and measured radius of a planet are 
wavelength dependent, and measuring the planetary radius 
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as a function of wavelength may allow the detection of ab- 
sorption features in the atmospheres, and thus enable obser- 
vations to infer the presence of atomic and molecular species 
(e.g. Seager fc Sasselov|2000 Brown|2001 1. Hot Jupiters are 
a class of planets with masses similar to Jupiter and on very 
short-period orbits, therefore they are intensely irradiated 
resulting in very hot atmospheres. Consequently, their large 
atmospheric scale heights make them particularly good tar- 
gets for this type of measurement. 

There are two approaches to measuring the wavelength 
dependence of a planet's radius. Transmission spectroscopy 
consists of monitoring the transit with spectroscopic mea- 
surements, which can be subsequently split into separate 
light curves for each wavelength channel, allowing a mea- 
surement of transit depth at each wavelength. This has 
been used to detect various species in the atmospheres of 



hot- Jupiters, including HD 209458 (e.g. Charbonneau et al. 
|2002} Na) and HD 189733 (e.g. |Swain et al.||2008[ H 2 OT 
CH4), and has also suggested a haze in the upper atmosphere 
of HD 189733 ( |Pont et al.|2008| ) at optical wavelengths. An- 
other method to measure the wavelength dependence of a 
planet's radius, is to make multi-colour photometric obser- 
vations of the transit light curve, and measure the transit 
depth for each. Spitzer observations of this type were used to 
infer the presence of H2O in the atmosphere of HD 189733 
( Tinetti et al.|200"7|, although this is disputed by Ehrenreich 



et al. (|2007|) and [Desert et al. (|2009 1 



These measurements have provided a wealth of infor- 
mation which feeds into theories of planetary atmospheres. 
However, both methods have advantages and disadvantages, 
and have been the subject of some controversy. For exam- 
ple, photometric measurements are much more straight for- 
ward, and the data reduction and analysis use established 
methods, and are therefore relatively robust against instru- 
mental systematic effects. However, observations are (typ- 
ically) taken during different transits, and stellar activity 
may cause variations of the measured planetary radius due 
to variable spots on the surface of the star, which change 
as the star rotates and reveals different hemispheres, and 
also evolve on timescales comparable to the period of typi- 
cal transiting planets ( jMosser et al.|2009| . 

Transmission spectroscopy on the other hand, avoids 
this problem by simultaneously monitoring the transit light 
curve at different wavelengths. Whilst there is still a chro- 
matic variance due to any spots on the stellar surface dur- 
ing a particular transit, this is a much smaller effect than if 
the surface spot distribution has changed. However, trans- 
mission spectroscopy is usually affected by systematic noise 
from the telescope and instrument (e.g. |Pont et al.||2008[ 
Swain et al. ||2008 1 , often larger than the signal we are try- 
ing to measure. This is commonly treated by constructing a 
multi-linear model of the baseline function, which describes 
how the measured light curve of the star would behave in 
the absence of a planetary transit due to changes in the 
state of the optics and detector. This is built from a linear 
combination of parameters that describe the optical state 
of the system, such as the position of the spectral trace on 
the detector, the temperature of the detector, and in the 
case of Hubble Space Telescope (HST) observations, the or- 
bital phase of the observatory. The baseline function is de- 
termined from the out-of-transit observations, and then used 
to decorrelate the light curve by projecting the function to 



the in-transit observations; however, the physical origins of 
these linear decorrelation models are poorly understood, and 
as we will see later, the choice of model has significant effects 
on the output transmission spectra. 

HST transmission spectroscopy with NICMOS (Near 
Infrared Camera and Multi-Object Spectrometer) has led 
to some of the most detailed studies of exoplanet systems, 
but it in particular suffers from these systematics effects (e.g. 



Swain et al. 2008; P ont et al. ||2009p , which are comparable 



to if not larger than the expected signal due to molecu- 
lar absorption. Of the NICMOS observations currently in 
the literature, one group have claimed the detection of sev- 
eral molecules in the atmospheres of various hot Jupiters 
(e.g. |Swain et al.|2008| HD 189733; [Tinetti et al.|20lo| XO- 
1). However, other analyses of NICMOS transmission spec- 
troscopy, have been unable to untangle the instrumental sys- 
tematic effects from wavelength dependent absorption from 
the planet (e.g. |Pont et al.||2009[ |Carter et al.||2009| >. Fur- 



thermore, for HD 189733, there have been significant dis- 
agreements with wavelength dependent photometric mea- 
surements of the light curve and transmission spectroscopy, 



with ISing et al. ( 2009 \ failing to detect water in the atmo- 



sphere of HD 189733 previously reported by | Swain et al.| 
( 2008 1 , using narrow band photometric measurements with 
NICMOS. [Tinetti et aL] ( |2010[ ) attribute this to variable spot 
distributions for each of transit, arguing multi-colour photo- 
metric observations cannot reach the required accuracy for 
molecular spectroscopy, despite |Sing et ah] ( |2009[ ) account- 
ing for this using continuous ground-based monitoring of the 
stellar flux. 

As these results are vital to our understanding of plan- 
etary atmospheres, we must explore the possibility that the 
wavelength dependance of the planetary radius is no more 
than systematic noise in the detector. As the HST orbits 
with a period of ~96 minutes, transit light curves are split 
into discrete orbits, typically sampling several out-of-transit 
portions of the light curve, used to establish the baseline 
function, and 1-2 in-transit portions. Trends are typically 
seen in each orbit, related to the orbital phase of the HST 
and its effects on the position of the spectral trace. Fur- 
thermore, as the optical system cannot be reset to exactly 
the same configuration for every orbit, systematic offsets be- 
tween the flux levels in each orbit may also result. Efforts 
to gain a deeper understanding of NICMOS systematics and 



obtain photon limited photometry arc ongoing (e.g. Burke 



et al. 2010). However, the offsets are not satisfactorily ad- 



dressed in much of the literature, as any systematic offset 
may not be corrected for or indeed may be over-corrected 
for when the baseline function is projected to the in-transit 
portion of the light curve in the case of an imperfect baseline 
model. These in-transit offsets will be hidden from the light 
curve residuals when fitting the transit depth, and therefore 
not taken into account in the error analysis. 

Fig. [l] shows a plot of the NICMOS spectra of transit- 
ing planets reported in the literature to dat^J transmission 
spectra for HD 189733 ( |Swain et al~1|2008|>, GJ-436 ( |Pont 



et al.||2009] ), and XO-1 ( jTinetti et ai.||2010[ ), and dayside 



1 For XO-1 and HD 209458 secondary, only every 5th point is 
plotted for comparison, because a boxcar smoothing was applied 
to the published data. 
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Figure 1. Plot of some of the NICMOS transmission and emis- 
sion spectra reported in the literature; transmission spectra for 
HD 18 9733 (|Swain et al.|2008| , GJ-436 | |Pont et a.l.]|2009| ) , and 
XO-1 ( |Tinetti et al.|201of and dayside emission spectra for HD 
189733 (Swain et al.|2009b| and HD 209458 (Swain et al.|2Q09a[ ). 
Each spectrum is plotted as % absorption as a function of wave- 
length, but offset for clarity. The peak-to-peak amplitudes of all 
reported spectra are approximately equal, and also show similar 
features. 



emission spectra for HD 1897 33 (Swain et al.|2009b| ) and HD 
209458 (Swain et al. 2009a I. The emission spectra are ob- 



tained in much the same way as the transmission spectra, by 
observing an eclipse of the planet, and using the same linear 
decorrelation methods to correct for the baseline function. 
It is clear that all five spectra show approximately the same 
peak-to-peak amplitude, and have wavelength dependance 
on similar scales. This is perhaps surprising, given the plan- 
ets have vastly different masses, radii and irradiation levels. 
Transmission and emission spectra also show the same am- 
plitudes. Whilst this is by no means evidence that the spec- 
tra are the result of systematic noise in the detector, and 
the similarity could be coincidental, it is nonetheless wise to 
explore this possibility. It is therefore important that the re- 
sults reported from NICMOS transmission spectroscopy be 
re-analysed to confirm previous detections of molecules, in 
particular because the results are dependent on a complex 
decorrelation procedure, with an arguably ad hoc model. 

Here, we present a re-analysis of the three NICMOS 
transmission spectra shown in Fig. [I] Sects. [2] [3] and [4] de- 
scribe the data reduction, analyses and results of the HD 
189733, GJ-436, and XO-1 data, respectively. Finally, Sect. [5] 
summarises and presents our conclusions. 



2.1 Observations 

A transit of HD 189733 was monitored on 25 May 2007 with 
HST/NICMOS, using the G206 grism covering the wave- 
length range 1.4 - 2.5 um. Analyses of these data were first 
reported in |Swain et aL] ( |2008| hereafter S08). As HD 189733 
is not in the continuous viewing zone of the HST, the tran- 
sit was observed over five half-orbits (~ 48 minutes each), 
consisting of 638 spectra in total, all with exposure times 
of 1.624 seconds. The first, second, fourth and fifth orbits 
cover the out-of-transit part of the light curve and consist 
of 119, 128, 131 and 130 observations, respectively. These 
data are required to determine the photometric baseline. 
Only the third orbit was taken in-transit, and consists of 
the remaining 130 images. In addition to the spectra, some 
exposures were taken at the beginning of the first orbit to 
enable wavelength calibration. 

The calibrated images, one of which is displayed in 
Fig. [2] include all basic calibrations except for flat-fielding. 
Correctly flat fielding spectroscopic data requires taking 
into account the wavelength dispersion and position of the 
source, and therefore a different flat-field correction would 
be required at each point along the spectral trace of the tar- 
get. S08 argue that flat-fielding is not required for differential 
spectroscopy, as each light curve will be normalised before 
measurement of parameters from the light curves. However, 
strictly speaking flat-fielding is required to accurately de- 
termine the background value in each wavelength channel. 
This is particularly important, as underestimating or overes- 
timating the background results in a variable transit depth, 
which is exactly what we are trying to measure. Further- 
more, the G206 grism has a relatively high background com- 
pared with NICMOS's other two grisms (about 800 electrons 
per pixel per exposure), due to thermal background emission 
from HST. Whilst the background is very stable temporally, 
unfortunately, it does significantly vary spatially over the 
detector by as much as 200 electrons per pixel per expo- 
sure. This will particularly effect the light curves extracted 
from the edge of the spectrum, where the overall counts 
are much lower. We did try flat-fielding each image using 
a flat-field taken with the G206 grism in place. This cor- 
rects for diffuse light over the whole detector. This will not 
effect the depth of the final normalised light curves except 
through the background value. This process does not pro- 
vide a satisfactory background correction, and the analyses 
described below were attempted with and without flat-field 
corrections. We also note that the background varies quite 
smoothly, and therefore is unlikely to be responsible for any 
narrow 'features' seen in the transmission spectrum, but can 
certainly affect the overall depth and shape. 

Prior to extracting the spectra, we identified the defec- 
tive pixels flagged by the calibration pipeline, and additional 
significant outliers, and corrected for them by replacing the 
pixel value with an interpolation of the surrounding eight 
pixels. A relatively small number of pixels required correc- 
tions, and in fact this process has very little effect on the 
output spectra. A ID spectrum was then extracted as follows 
for each image. For each column along the spectral trace in 
the dispersion axis (x), a centroid along the spatial axis (y) 
was calculated to determine the position of the spectrum. A 
sum of 35 pixels along each column centred on this position 
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Figure 3. Extracted ID spectra of a typical in-transit (green) 
and out-of-transit (blue) observation of HD 189733 showing the 
number of electrons collected per pixel channel. The wavelength 
decreases with increasing x position. 



Figure 2. NICMOS image of HD 189733 taken with the G206 
grism. The blue line marks the centroids of the first-order spec- 
trum of HD 189733, with the green lines flanking representing the 
extent of the extraction region. Immediately left is the second- 
order spectrum, and below is a companion star to HD 189733. 
The bright band at the right of the detector is a feature of G206 
grism spectra, caused by the warm edge of the aperture mask. 



was used to extract the flux for each wavelength channel, 
after subtraction of a background value for each pixel. A 
width of 35 pixels was chosen to minimise the RMS in the 
white light curve of orbits 2, 4 and 5. This was repeated for 
90 pixel columns along the dispersion axis. The extraction 
regions used are marked in Fig. [2] For the background sub- 
traction, we experimented with various techniques. First, a 
global background subtraction was used, similarly to S08. 
The background was taken to be the average of a large un- 
illuminated region above the spectral trace (we tested using 
different regions). This is not ideal, as previously mentioned 
the background varies spatially over the detector. As a first 
order correction, we instead calculated the background sep- 
arately for each pixel column, as the average value of the 
un-illuminated region above the spectral trace along the col- 
umn. Again we note that this does not provide a satisfac- 
tory correction, as the background varies along both x and y. 
Both global and wavelength dependent corrections were used 
for the subsequent analysis. For the remainder of this paper, 
we will display results from data that is not flat-fielded, and 
using separate columns for background subtraction, but our 
conclusions remain the same for each case. 

Extracted spectra from a typical in-transit and out-of- 
transit observation are shown in Fig.[3j giving approximately 
430 000 electrons in the brightest pixel channel, and approx- 
imately 120 000 electrons in the faintest channel. The fea- 
tures in the spectra do not correspond to real stellar fea- 
tures, but result in variation of the sensitivity of the grism 
with wavelength. Each wavelength channel in the ID spec- 
tra is then used to construct a time-series, after binning in 
5 pixels along the dispersion direction, resulting in 18 light 
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Figure 4. Raw 'white' light curve of HD 189733 obtained by in- 
tegrating the flux from each spectra over all wavelengths, showing 
the sampling of the transit, alongside its best-fit model. While all 
orbits are seen to suffer from systematics, orbit one exhibits by 
far the largest, commonly attributed to spacecraft 'settling', and 
is excluded from subsequent analysis. 



curves. A 'white' light curve was constructed by integrating 
the flux over the entire wavelength range for each image. 
This was extended to include 110 pixel columns, which min- 
imises systematics arising from small changes in the position 
of the spectral trace. This is plotted in Fig. [4] and shows the 
sampling of the light curve over the five orbits. Systemat- 
ics are evident in each orbit, but particularly for orbit one. 
This is commonly found in similar NICMOS data (see e.g. 



Sects. 3.1 and 4.1|. It is attributed to spacecraft 'settling' 
(e.g. S08; |Pont et al.|2009[ ), and orbit one is excluded for the 
remainder of this work. 

The raw light curves are shown in Fig. [5] for each of the 
18 wavelength channels, after normalising by fitting a line 
through orbits 2, 4 and 5, and dividing the light curves by 
this function. The light curves are clearly seen to exhibit 
strong time-correlated noise, which needs to be removed if 



© 2002 RAS, MNRAS 000,[Tp7l 



A new look at NICMOS transmission spectroscopy 5 



l *H*f*»* t - - 











lliWKtfJK 1 - - - 


V 




, - - 


— mtHi** 








/ 


— -r*i)ri»*ai»- 


- - 


/ 


, -**jf*4**»* - - 








/ 


tH'.t*Ji*i 


lf«Wlt**»- , 




,- -«tf*#**»*- - 






/ 


/ 




jMtHHt* ■ 






r -tkft*4t*f*- - - 


- - -fft.fHt***- - — 




. ru!+n*$**- - 




/ 




t»4*f*.H« 


/ 


, »*- - - 


*<***1t*1* 





. **>***«?«» - 










/ 




********** 






/ 






**"W!,W(M 






/■ -4»t«»«*ir — 


*<J*ri*i-»r 




. »Wfirt»« - 








- — -lif**i\**i+- - - 


/ 


, - 


-</*«VI«*»» 

















, *t*f*t**f - - 


********** 


»i*tMV»»«»- 






s - - 


»rtJ***W» 


«I4*Kl„ fJ - 






, JUjU Htr - 


***'**'T»* 










««H«**»M 


»;I*ik*Vt#f< 




/■ - - 








, -*MW«»f - - 




4*it«4fMt* 




ewttrf**?* 


- ^H.o,r - 






, J«iMf»ir - 


- - 


- - -frfc«>H«r - - 


^ - 




, _4KW«^ fr _ 






- - -f<*ft**^> 


- - -#*M**(hsi 








#s^ii*4 S( 


- - -j^J*^>r - - 



















0.02 
0.00 
-0.02 
-0.04 
-0.06 
-0.08 



65.12 
2 65.08 
~ 65.04 
< 65.00 
64.96 
2.92 
3 2.88 

Q. 

g 2.84 

2.80 

2.00 
1.98 
1.96 
1.94 
f 1.92 
1.90 
1.88 
77.16 



246.05 246.10 246.15 246.20 246.25 
[HJD] - 2454000 



246.30 



77.14 
g 77.12 ■ 
H 77.10 

77.08; 

246.05 



/ / / / 



246.10 246.15 246.20 246.25 
[HJD] - 2454000 



246.30 



Figure 5. Raw light curves of HD 189733, for each of the 18 
wavelength channels from 2.50 u.m (top) to 1.48 urn (bottom), af- 
ter normalising by fitting a linear function through orbits 2, 4 and 
5. Orbit 1 is not plotted as it is discarded for the analysis. Clearly, 
time-correlated noise is present in each light curve, and must be 
removed to measure the transmission spectrum. The dashed grey 
lines show a transit model generated for HD 189733, used to guide 
the eye. 



we are to measure the transmission spectrum at the level re- 
quired to detect molecular features. The systematics may be 
understood to arise from small motions of the spectra across 
the detector, related to the orbital motion of the HST. Refer- 
ring back to Fig. [3j the wavelike features in the spectra will 
move into different wavelength channels and cause features 
on short wavelength intervals, even if the stellar spectrum 
and flat-fielding are smooth. 

The optical state parameters, described by S08, were 
therefore measured in an effort to model and remove these 
systematic effects. We extracted the shift of the spectral 
trace along the x axis (AX) by cross-correlation of the ID 
spectra, the shift of the position of the y axis by averaging 
the centroids determined earlier (AY), and the angle the 
spectral trace makes with respect to the x-axis (8) by fitting 
a line to the centroids of the spectra. The width (W) of each 
spectrum was also measured by fitting gaussian functions 
along each extraction column, and taking the average. In 
the absence of a direct measurement of the temperature (T) 
of the detector, or a proxy calculated from the bias levels 
( Gilliland fc Arribas|2003| ), a proxy for this was taken as the 
temperature of the NIC1 mounting cup (S08). However, this 
is not monitored at a high enough precision to accurately 
monitor the temperature. The temperature is an important 



Figure 6. Extracted state parameters for HD 189733 observa- 
tions for orbits 2-5. For orbit 3 we reconstruct the model baseline 
as a linear function of all state parameters, plus a constant flux 
level, and the orbital phase of the HST and its square. The co- 
efficients for each parameter of the model are determined from 
orbits 2, 4 and 5. 



parameter to describe the state of a NIR detector, and this is 
likely an important limitation of NICMOS analysis. A plot 
of the optical state parameters and the detector temperature 
is shown in Fig. [6] 

It is important to note that we extract our parameters 
slightly differently from S08 in some cases. In particular we 
used an averaged value for the shifts in position AX and 
Ay per image, in order to monitor the overall optical state 
of the detector, and therefore had identical parameters for 
all wavelength channels, unlike S08. S08 determined x and 
y positions by fitting spectra along diagonals and rotating 
back along the two nominal axis, resulting in different pa- 
rameters for each wavelength channel, and in some cases 
different trends. In either case, the parameters show a sim- 
ilar dispersion and amplitude, and should allow a similar 
removal of the systematics. 



2.2 Analysis 

Each light curve was decorrelated using the multi-linear 
baseline model from S08. This assumes that the baseline 
function can be constructed by a linear combination of the 
optical state vectors AX, AY, W and 9, the proxy for the 
detector temperature T, and also the orbital phase of the 
HST (4> H ) and its square (<^|/). In addition, we used a con- 
stant vector to represent the out-of-transit flux level (/o), 
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as clearly the overall flux level is not given by a linear func- 
tion of the state variables, but rather the state variables 
are used to correct for variations from the expected out- 
of-transit flux. These parameters are hereafter collectively 
referred to as the 'decorrelation parameters'. For each wave- 
length channel, the out-of-transit baseline flux is given by 

N 

Vi = /J PkXi t k + e< 
fe=i 

at time i, where /?& are the coefficients for each of the N 
decorrelation parameters (in this case; fo, AX, Ay, W, 8, 
T, 4>h and (f>%), X it ^ is the value of each decorrelation pa- 
rameter k for time z, and e, the corresponding residual. This 
may be easily written in matrix form as 

y = X/3 + ?, 

where X is the state matrix containing the measured de- 
correlation parameters, f3 is a vector containing the coeffi- 
cients of each decorrelation parameter, and e are the residu- 
als. The best-fit coefficients j3 are then found by linear least 
squares 

$ = (X T X)- 1 X T y, 

for orbits 2, 4 and 5 only. The model baseline function b can 
then be reconstructed for all orbits from a linear combina- 
tion of P and X; 

b = X/3. 

Each light curve is then decorrelated by dividing through 
by the model baseline function. This technique can only be 
expected to work if we interpolate for the in-transit orbit, 
and if the baseline function is well represented by the linear 
model over the range of the decorrelation parameters. 

For this analysis, each light curve is treated indepen- 
dently, with j3 calculated separately for each. An example of 
this decorrelation process on one of the wavelength channels 
is shown in Fig. [7] and all of the decorrelated light curves 
are shown in Fig. U S08, who use a similar procedure to 
this, conclude that it satisfactorily removes time-correlated 
systematics from the light curves. Whilst our results agree 
with this conclusion for the out-of-transit orbits, residual 
systematics are clearly visible in the in-transit orbit, which 
is the important one for determining the radius ratio. 

In order to obtain the transmission spectrum from the 
light curves, a transit model must be fitted to each decorre- 
lated light curve to measure the planetary radius. We used 
the transit model described in | Gibson et al.| ( |2008[ ), which 
assumes a circular orbit to calculate the normalised sepa- 
ration of the star and planet centres (z) as a function of 
time, using the orbital parameters and the masses and radii 
of the star and planet. The analytic models of Ma ndel "&] 
Agol ( 2002 1 were then used to calculate the light curves 



from z, the ratio of the planet-to-star radii p, and the limb 
darkening parameters. The limb darkening parameters were 
calculated in each of the wavelength channels for a quadratic 
limb darkening law (D. Sing, private communication) using 
the methods described in Sing (20101. The orbital parame- 
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Figure 7. Example of the decorrelation procedure on one of the 
wavelength channels. The top plot is of the raw light curve in elec- 
tron counts. The green points represent the reconstructed baseline 
function (b) used to decorrelate the light curve. The bottom plot 
shows the decorrelated light curve along with its best-fit transit 
model, used to determine the planet-to-star radius ratio. Whilst 
the out-of-transit orbits are usually whitened, residual system- 
atic noise is often seen in the in-transit orbit, showing the linear 
baseline model is not robust. 



decorrelated light curve was then fitted for p, to obtain the 
transmission spectrum. 

To calculate the uncertainty in p for each wavelength 
channel, we used a residual permutation (or 'prayer bead') 
algorithm (see e.g. Gillon et al. 2007 Southworth 2008 1, sim- 



ilar to that used in Gibson et al. (2009 20101. This method 



ters of the system and the stellar mass and radius were held 
fixed at the values given by |Pont et al.| ( [2008] ), and the cen- 
tral transit time was determined from the ephemeris. Each 



accounts for the correlated noise in the light curve by re- 
constructing light curves by combining the best-fit model 
and its residuals, but shifting the residuals before combin- 
ing them to determine the effects of the correlated noise on 
the determined parameters. This method preserves both the 
correlated and random noise in the resampled light curves, 
and therefore both are taken into account when determining 
uncertainties. 

The residual permutation was applied to the raw light 
curves prior to performing the decorrelation procedure, in 
order to take into account uncertainties from both the linear 
decorrelation and the light curve fitting. The best fit baseline 
function b and transit model rh^f were determined as before, 
and the residuals from the fit f are used to reconstruct the 
light curve, but each time the residuals are shifted by a 
random offset to give f p . Any residuals that fall off the 'edge' 
are looped back to the beginning. The new light curve y p is 
then reconstructed by adding the shifted residuals to the 
best-fit model, followed by a pointwise multiplication of the 
best-fit baseline function; 

Hp = ("W + r p ) ■ (X/3) 

The decorrelation and light curve fitting is done as before, 
to determine p. The procedure is repeated 1000 times with 
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Figure 8. Decorrelated light curves for each of the 18 wavelength 
channels. Much of the time-correlated systematic noise is removed 
from the light curves, as seen from a direct comparison with Fig. [5] 
However, some of the correlated noise remains particularly in the 
in-transit orbit, and must be taken into account in the error anal- 
ysis of the transmission spectrum. The dashed grey lines again 
show a transit model generated for HD 189733. 



random perturbations to the shift in residuals, each time 
varying the starting value for p to ensure the starting pa- 
rameters do not affect the results. The resulting distribution 
of p is then used to estimate its uncertainty in each wave- 
length channel. 



2.3 Results 

The resulting NICMOS transmission spectrum for HD 
189733 is shown in Fig. [9] The data from S08 are also plot- 
ted for comparison, after converting from transit depth to p. 
For the most part, the spectra show the same basic shape. 
It is not clear exactly where the discrepancies in a few of the 
wavelength channels arise, but they are likely explained by 
one or a combination of the following; different pixel columns 
and widths used for the wavelength channels, different meth- 
ods used to determine the background, the fact that we fit 
for the light curve rather than just taking an average of 
the in-transit orbit, that the decorrelation parameters are 
extracted slightly differently, and finally the corrections ap- 
plied by S08 for limb darkening and star spots. We were 
unable to reproduce exactly the same results as S08 using 
a global background correction. Another obvious difference 
between the two spectra is that the uncertainties we calcu- 
late using the residual permutation method are significantly 
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Figure 9. Transmission spectrum of HD 189733 generated by 
determining the planet to star radius ratio from the decorrelated 
light curves for the 18 wavelengths channels. Our results are the 
red points, and those from S08 are shown in grey for comparison. 
The horizontal dashed line is the radius ratio measured for the 
white light curve. 



larger than those given in S08, particularly at the edge of 
the spectrum. 

However, we do not believe the residual permutation 
method fully accounts for the uncertainties in p. Fig. [10] 
shows a plot of the in-transit residuals for three wavelength 
channels. The residuals are of decorrelated light curves, af- 
ter subtracting a model generated from the planet-to-star 
radius ratio measured from the white light curve, and us- 
ing limb darkening parameters specific to each wavelength 
channel. The difference between models generated for the 
best-fit planet-to-star radius ratio for each channel (i.e. the 
value in the transmission spectrum), and the planet-to-star 
radius ratio from the white light curve (i.e. a constant ra- 
dius model) are also shown. As is clear for all three channels, 
systematic noise is present at a level comparable to or larger 
than the difference between the two models. The deviations 
from a constant transmission spectrum may therefore arise 
from systematics not removed from the in-transit orbit. A 
similar level of systematics is visible in the residuals for all 
wavelength channels. 

We are yet to address one final correction on each wave- 
length channel carried out by S08, the 'channel-to-channel' 
corrections. This involves taking the weighted average of the 
light curve residuals for all the wavelength channels, and 
subtracting them from individual channels after decorrelat- 
ing the light curves. This should remove any common time- 
correlated systematics from the light curves. To check that 
this did not affect our results significantly, we carried out 
this step. Fig. |11| shows the same light curves as Fig. [To] af- 
ter the channel-to-channel correction has been applied. This 
process does marginally reduce the RMS of the residuals, 
but does not completely remove the systematics. This is not 
surprising given all wavelength channels do not show the 
same structure of systematics. We would also not expect 
this to change the transmission spectrum, as applying the 
same correction to all light curves will only shift the depths 
uniformly. However, reducing the RMS may have an impact 
on the uncertainties. The resulting transmission spectrum 
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Figure 10. Residuals of the in-transit orbit from a 'constant ra- 
dius' model generated from the white light curve planet-to-star 
radius ratio for three wavelength channels; 2.41, 2.07 and 1.66 
urn from top to bottom, respectively. The limb darkening coeffi- 
cients used were specific to each wavelength channel. The hori- 
zontal black dashed line shows the constant radius model, and the 
dashed red line the difference between it and the best-fit model 
determined when fitting p for each wavelength channel. It is clear 
that systematic noise is present in the in-transit orbit for the three 
channels, at a level comparable to or even larger than the devia- 
tions from the constant radius. A similar level of systematics can 
be seen in all wavelength channels. 



is shown at the top of Fig. |12| This does not cause signif- 
icant changes in the transmission spectrum, or even in the 
calculated uncertainties. This indicates that the uncertain- 
ties in the transmission spectrum are probably dominated 
by determination of the baseline function. Recently, |Burke| 
et al. (20101 identified systematic effects that should be ac- 



counted for with NICMOS grism spectroscopy, in particular 
'gain-like' variations that arise from seven states of the de- 
tector electronics. Importantly, they noted that the channel- 
to-channel correction would likely account for this. 

Significant uncertainties in the depth will also come 
from offsets in the baseline flux level between the orbits. Off- 
sets in the flux levels are seen in the raw light curves, and 
could easily be induced by the decorrelation model (as well 
as corrected for), as the decorrelation parameters themselves 
are seen to have offsets between orbits, therefore clearly a 
linear baseline model can produce 'artificial' offsets. These 



flux offsets were also identified by Carter et al. (20091, who 



tried to fit decorrelation models based on state parameters 
to attempt to remove the flux offsets, and concluded the 
procedure was unreliable with no physical justification for a 



Figure 11. Same as Fig. |10| but after the channel-to-channel cor- 
rection has been applied. In some cases the RMS of the residuals 
is reduced, but there are still systematics present that are compa- 
rable to or even larger than the deviations from a constant radius 
model. 



linear baseline model. As variations in the flux level for the 
in-transit orbit are fitted for when determining the plane- 
tary radius, the residual permutation algorithm cannot fully 
take this additional uncertainty into account, as the resid- 
uals will not contain the signal of any induced offset. To 
test for the possibility that the decorrelation methods may 
not fully correct for the flux level, we repeated the decor- 
relation process only using orbits 2 and 4 to determine the 
decorrelation coefficients. Referring to Fig.|6j we are still in- 
terpolating for the in-transit orbit for all the decorrelation 
parameters. Therefore, if the model for the baseline flux is 
correct, we should still see similar results. Fig. |13| shows an 
example of this decorrelation process. It is clear that orbit 5 
has a significant offset from the baseline flux of orbits 2 and 
4. This shows that the linear decorrelation method is not ro- 
bust when extrapolating over the decorrelation parameters. 
If the functions describing the baseline function are not lin- 
ear in nature (which is by no means clear, as even assuming 
the baseline flux can realistically be written in terms of the 
optical state parameters is a strong assumption), a spurious 
offset could just as easily be induced in the in-transit orbit, 
leading to channel dependent transit depths, and resulting 
in a variable transmission spectrum. 

The transmission spectrum obtained when decorrelat- 
ing only using orbits 2 and 4 is shown in Fig. |12| As expected, 
the uncertainties are now significantly larger than before (or- 
bit 5 is not used in the fitting process). However, the trans- 
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Figure 12. Transmission spectra of HD 189733 generated by de- 
termining the planet radius p for the 18 wavelengths channels. 
The top plot shows the transmission spectrum obtained after ap- 
plying the channel-to-channel correction. The second plot shows 
the spectra by decorrelating the light curves only using orbits 2 
and 4. The third plot shows the result of a 'quadratic' decorrela- 
tion, and the bottom plot shows a linear decorrelation but with 
the angle vector (8) removed from the state matrix (X). Again, 
our results are the red points, and those from S08 are shown in 
grey for comparison. The horizontal dashed grey line marks the 
level of the planet-to-star radius ratio of the white light curve. 
For the last three cases the spectra are altered considerably from 
the first, and none of the features reported in S08 are present in 
all spectra using the different methods. 



mission spectrum shows a significantly different structure 
than before (c.f. FigS|. In fact the only 'feature' common to 
both, is the dip at around 2.1 (j.m, although its depth relative 
to the uncertainties varies. The edges of the spectrum show 
quite different behaviour. We interpret the significant differ- 
ences in the two spectra as evidence that the linear baseline 
model is not sufficient to correct the light curves, and can 
induce unwanted offsets in the in-transit orbit, which are 
indistinguishable from a real atmospheric signal. 
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Figure 13. Example of the decorrelation procedure on one of the 
HD 189733 wavelength channels, similar to Fig. [7] but showing 
a different wavelength channel, and now decorrelating only using 
orbits 2 and 4. It is clear that the decorrelation does not prop- 
erly correct for orbit 5, and an artificial offset is introduced by 
the linear baseline model. Note that this offset is larger than the 
features in the reported transmission spectrum. We argue that 
the decorrelation process may introduce similar spurious offsets 
to the in-transit orbit, indistinguishable from a real atmospheric 
signature. 



We also conducted additional tests to see how de- 
pendent the transmission spectrum is on the decorrelation 
model. A 'quadratic' decorrelation was tried by adding the 
terms AX 2 , AY 2 , W 2 , 8 2 and T 2 , to the state matrix (<j> H 
already has a quadratic term). We also made tests decorre- 
lating the spectra after removing each vector from the state 
matrix. We found that the most important decorrelation 
parameter to be 9, and to a lesser extent AY. The other 
decorrelation parameters play a much less important role. 
This is consistent with the findings in S08, and is rather 
obvious after closer inspection of the decorrelation parame- 
ters in Fig. [6] and AY show significant offsets between the 
in-transit orbit and the out-of-transit orbits. Thus, decorre- 
lating with these state parameters causes larger shifts in the 
in-transit orbit than the others, and consequently the most 
significant corrections to the transmission spectrum. 

The spectra produced using the 'quadratic' decorrela- 
tion, and after removing 9 from the state matrix are shown 
in Fig. |12| Note all of the spectra plotted in Fig. [12] are af- 
ter the channel-to-channel correction has been applied, but 
its exclusion does not alter the results significantly for any 
case. It is clear from these plots that the various techniques 
used to decorrelate the light curves produce quite different 
output transmission spectra. Each of the features reported 
in S08 vanishes in at least one of the cases. As 9 was pre- 
viously shown to be the most important parameter in the 
decorrelation, the interpretation of molecular features in the 
spectrum is heavily dependant on a linear dependance of the 
flux level with the angle the spectrum makes on the detec- 
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tor. We do not argue that any of these methods are better 
than the others or the method used by S08, but they give 
incompatible results and there is no way to distinguish be- 
tween them. Whilst most of the systematic effects in the 
light curves are caused by intra-pixel sensitivity changes as 
the spectral trace drifts and rotates on the detector, it is not 
clear that a linear decorrelation will remove these systemat- 
ics. In fact, we have shown that the decorrelation can cause 
spurious offsets in the in-transit flux level. 

Furthermore, the largest disagreements between the 
spectra occur at the edges, where the background is expected 
to significantly affect the transit depth. The ratio of the true 
transit depth (D) is related to the measured transit depth 
(D') by 

D 



1 + H' 



where B is the (uncorrected) background counts, and Fq the 
out-of-transit counts. Therefore, an underestimated or over- 
estimated background will lead to incorrect transit depths. 
This effect scales with the ratio of the background to the 
out-of-transit flux. Thus where the flux count is lower, i.e. 
at the fainter channels (the edges) of the spectrum, this ef- 
fect will be much greater. For the faintest channels, we have 
approximately 150 000 x 5 electrons (Fig. [3| collected over 
35 x 5 pixels. If we assume a background correction error of 
50 electrons per pixel (a conservative estimate considering 
the background can vary by up to 200 electrons per pixel), 
the ratio of true depth to measured depth is about 1.012, 
more than enough to significantly affect the measured tran- 
sit depth. Clearly we cannot trust the edges of the spectra 
to give us an accurate measurement of the planet-to-star 
radius ratio. 



3 GJ-436 

3.1 Observations 

A transit of GJ-436 was monitored on 10- 11th November 
2007 with HST/NICMOS, using the G141 grism covering 
the wavelength range 1.1 - 1.9 um. Analyses of these data 



were first reported in Pont et al. (2009 hereafter P09). We 
only used the first of two light curves reported in P09 for 
this work, as they concluded that the noise and systematic 
effects were lower during the first visit. GJ-436 is again not 
in the continuous viewing zone of the HST, and the transit 
was observed over four half-orbits, consisting of 935 spec- 
tra in total, all with exposure times of 1.993 seconds. The 
first, second and fourth orbits cover the out-of-transit part 
of the light curve and consist of 180, 252 and 251 observa- 
tions, respectively, and are used to determine the photomet- 
ric baseline. The third orbit was taken in-transit, covering 
mid-transit and egress, and consists of the remaining 252 
images. Similarly to HD 189733, some exposures were taken 
at the beginning of the first orbit to set the wavelength cal- 
ibration. 

Again, calibrated images include all basic calibrations 
except for flat-fielding. The background with the G141 grism 
is much lower than with the G206 grism, and is therefore the 
background correction is less of an issue for this data set. The 
background was again very stable over the duration of the 
observations, with a typical level of only 25 — 30 electrons 
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Figure 14. Extracted ID spectra for GJ-436 of a typical in- 
transit (green) and out-of-transit (blue) observation showing the 
number of electrons collected per pixel channel. The wavelength 
decreases along the x-axis. 



per pixel. The spatial variation over the detector is about 
10 electrons per pixel. We again tried flat-fielding using a 
flat-field taken with the G141 grism in place, but this did 
not provide a satisfactory correction and we did not use the 
flat-field corrected images for this analysis. 

The spectra and light curves were extracted using the 
same technique as for HD 189733. The spectra were ex- 
tracted along 90 pixel columns, with a width of 20 pix- 
els. The width of 20 pixels was selected to minimise the 
RMS in orbits 2 and 4 of the white light curve. Again we 
experimented using both global and column-specific back- 
ground corrections, but it made insignificant difference due 
to the low background. Fig |14| shows a ID spectrum for 
an in-transit and out-of-transit observation. About 420 000 
electrons were collected per exposure in the brightest pixel 
channel, and approximately 170 000 electrons in the faintest 
pixel channel. The raw light curves were then constructed as 
for HD 189733 by binning the spectra into 5 pixel bins, and 
are shown in Fig. [15] They show similar systematic effects 
to HD 189733. The white light curve was extracted using 
110 pixel columns and is shown in Fig. |16| excluding orbit 1 
(left out of the remaining analysis due to large systematics), 
and the optical state parameters are shown in Fig. [17] 



3.2 Analysis 

Each light curve was decorrelated using the procedure de- 
scribed in Sect. [2~2l The decorrelation coefficients were de- 
termined for orbits 2 and 4. Referring to Fig. |17| it is clear 
that we are extrapolating for the angle of the spectrum when 
correcting the in-transit orbit. We therefore must exclude 
this parameter from the multi-linear decorrelation, other- 
wise this results in a spectrum with very large features 
caused by spurious offsets in the in-transit baseline func- 
tion. An example of the decorrelation procedure is shown in 
Fig. [18] 

The decorrelated light curves are shown in Fig. |19| Each 
of the light curves was then fitted using the transit models 
described in Sect. |2.2] to determine p as a function of wave- 
length, fixing the stellar and orbital parameters at those de- 
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Figure 15. Raw light curves for each of the 18 wavelength chan- 
nels for GJ-436, from 1.87 Jim (top) to 1.19 |J.m (bottom), after 
normalising each one by fitting a linear function through orbits 
2 and 4. Similarly to HD 189733, correlated noise is seen in each 
light curve. 



termined by P09, who undertook a thorough analysis of the 
integrated light curve. The central transit time was deter- 
mined from the white light curve (as a well sampled egress 
allows an accurate determination of the central transit time) , 
and fixed at this value when fitting for each wavelength chan- 
nel. The limb darkening parameters were calculated using 



the method of Sing (20101 as for HD 189733, and the uncer- 



tainties were calculated using the residual permutation algo- 
rithm described previously, with 1000 light curves generated 
with random shifts applied to the light curve residuals. 



3.3 Results 

The resulting transmission spectrum for GJ-436 is shown 
in Fig. |20| There seems to be a significant variation of 
the radius ratio with wavelength. Unfortunately, with only 
two out-of-transit orbits to work with, we do not have 
enough data to test the decorrelation model with differ- 
ent orbits. The planet-to-star radius ratio is calculated as 
p = 0.0830 ± 0.0005 for the white light curveQ consistent 
with P09. The transmission spectrum does not appear con- 
sistent with this value, as nearly all wavelength channels 
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Figure 16. Raw 'white' light curve of GJ-436 found by integrat- 
ing the flux from each spectra over all wavelengths, showing the 
sampling of the GJ-436 transit, excluding the first orbit. 
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Figure 17. Extracted de-correlation parameters for GJ-436 dur- 
ing orbits 2—4 as a function of time. For the angle of the spectral 
trace 8, the baseline function must be extrapolated for the in- 
transit orbit, hence this parameter cannot be used in the decor- 
relation process, without introducing large offsets in the transit 
light curve. 



2 This uncertainty is obtained when only fitting for the transit 
depth and a linear baseline function, and therefore is likely un- 
derestimated. 
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Figure 18. Example of the decorrelation procedure on one of the 
GJ-436 wavelength channels, showing the raw and decorrelated 
light curve at the top and bottom, respectively. The green points 
represent the baseline function. 
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Figure 19. De-correlated transit light curves of GJ-436 for each 
of the 18 wavelength channels. Again, much of the time-correlated 
systematic noise is removed from the light curves, but some re- 
mains. 
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Figure 20. Transmission spectrum of GJ-436 generated by deter- 
mining the planet-to-star radius ratio from the decorrelated light 
curves for the 18 wavelengths channels. The horizontal dashed line 
shows the planet-to-star radius ratio calculated from the white 
light curve. 



in the transmission spectrum give a larger depth. However, 
the transmission spectrum is consistent with the white light 
curve if we decorrelate only using the orbital phase parame- 
ters. A likely explanation for this is that the linear decorre- 
lation process induces false offsets in the in-transit flux level, 
as was suspected for HD 189733. 

It is interesting that the amplitude of the measured 
transit depth (~ p 2 ) shows a similar amplitude to that of 
HD 189733. Any variations in the transmission spectrum is 
likely due to systematic noise in the detector that we have 
not accounted for in our error analysis, rather than any real 
physical effect. Perhaps this represents the limit of NIC- 
MOS transmission spectroscopy observations due to corre- 
lated noise, without a more robust way to remove or control 
the systematics. 



4 XO-1 

4.1 Observations 

A transit of XO-1 was observed on 21 February 2008 with 
HST/NIC MOS, using the G141 g rism. These data were first 
reported in |Tinetti et al.| ( |2010] hereafter T10). XO-1 is a 
similar spectral type to HD 189733, but it is considerably 
fainter, with a J-band magnitude of 9.94 (c.f. 6.07 for HD 
189733). The transit was observed over five half orbits, con- 
sisting of 279 spectra in total, each with an exposure time of 
40 seconds. The first, second and fifth orbit cover the out- 
of-transit part of the light curve, and consist of 56, 56 and 
55 spectra, respectively. The third orbit covers the ingress 
with 56 spectra, and the fourth covers mid-transit, consist- 
ing of 55 spectra. XO-1 has a longer transit duration than 
HD 189733 or GJ-436, hence two consecutive orbits cover 
in-transit data. Again, some exposures were taken prior to 
the first orbit to set the position for wavelength calibration. 

The images were treated and the light curves extracted 
using the techniques described in Sects. [272] and [3~2| The ID 
spectra were extracted along 90 pixel columns of width 16 



© 2002 RAS, MNRAS 000,[TJ{T7] 



A new look at NICMOS transmission spectroscopy 13 




60 80 
x position (pixels) 



120 



Figure 21. Extracted ID spectra for XO-1 of a typical in-transit 
(green) and out-of-transit (blue) observation showing the number 
of electrons collected per pixel channel. 
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Figure 22. Raw 'white' light curve of XO-1 found by integrat- 
ing the flux from each spectra over all wavelengths, showing the 
sampling of the XO-1 transit, excluding the first orbit. 



pixels. The background was estimated from a column strip 
above each wavelength channel. We also experimented with 
flat-fielding and global background corrections, but neither 
significantly affects the final results. Fig|21| shows typical 
ID spectra extracted for in-transit and out-of-transit obser- 
vations, with the number of electrons collected per image 
per pixel column ranging from approximately 160 000 to 
400 000. Fig. |23| shows the light curves extracted for each of 
the 18 wavelength channels, after binning into 5 pixel bins. 
Fig. |22| shows the resulting white light curve, extracted over 
110 pixel columns. The first few points in each orbit exhibit 
strong variation in flux, and are thus excluded from the sub- 
sequent analysis. 

The white light curve shows much larger systematics 
than either the HD 189733 and GJ-436 light curves. Conse- 
quently, the light curves for each wavelength channel suffer 
from significant levels of correlated noise. The decorrelation 
parameters were extracted as for HD 189733 and GJ-436, 
and are plotted in Fig. |24| 
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Figure 23. Raw light curves of XO-1 for each of the 18 wave- 
length channels from 1.87 |im (top) to 1.19 (im (bottom) excluding 
orbit 1 , and after normalising each one by fitting a linear function 
through orbits 2 and 5. 



4.2 Analysis 

Given the large amount of systematics, one would like to 
apply the decorrelation used for HD 189733 and GJ-436 on 
these data. However, it is clear from Fig. |24| that for the 
in-transit orbits (3 and 4), we must extrapolate if we are 
to use corrections for AX, AY and 9. We therefore must 
exclude these decorrelation parameters from the procedure. 
An example of the decorrelation procedure using only W, T, 
b 2 H , is shown in Fig. 



and ( 



25 



The large residuals show the 



decorrelation process does not provide a more satisfactory 
correction. 

The decorrelated light curves are shown in Fig. |26| 
Time-correlated noise still remains in the light curves, but no 
combination of decorrelation parameters produces a satisfac- 
tory result. The light curves were then fitted using the tran- 
sit models described previously, fixing the orbital and stel- 



lar parameters at those determined by Torres et al. (2008) 



The central transit time was again set by fitting the white 
light curve with a transit model, and was then fixed when 
fitting the light curve for each wavelength channel. The un- 
certainties were determined using the residual permutation 
algorithm described, with 1000 realisations. 



4.3 Results 

Fig. |27| shows the resulting transmission spectrum for XO-1. 
The results from T09 are also plotted for comparison. The 
amplitude in terms of flux is >0.3 %, and the corresponding 
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Figure 24. Extracted optical state parameters for XO-1 during 
orbits 2-5. For the in-transit orbits (3 and 4), we have to ex- 
trapolate for AX, AY and 9, making the decorrelation process 
extremely difficult. 



Figure 26. De-correlated light curves of XO-1 for each of the 
18 wavelength channels. Some of the time-correlated systematic 
noise is removed from the light curves; however, a significant 
amount remains which could not be removed in the decorrelation 
procedure with any linear combination of decorrelation parame- 
ters. 
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Figure 25. Example of the decorrelation procedure on one of 
the XO-1 wavelength channels, showing the raw and decorrelated 
light curve at the top and bottom, repectively. The green points 
represent the baseline function. It is clear this does not provide 
a satisfactory correction for the systematics, which is obvious in 
the residuals for all orbits. 



uncertainties are about 0.05%. Both of these are significantly 
larger than those reported in T10, and obviously we have 
not reached the precision required to detect any molecular 
species. The reason for this is clearly related to the decor- 
relation parameters. Unfortunately, T10 did not provide a 
plot of their decorrelation parameters, or plots of the light 
curves. This makes it difficult to check where our reductions 
differ. As T10 did not mention the need to extrapolate for 
any of the decorrelation parameters, we must presume the 
difference between our analysis and theirs resides chiefly in 
the extraction of the decorrelation parameters. 



However, recently Burke et al. (20101 reported a broad- 



band analysis of this light curve, and provided plots of 
AX, AY and 6. The slight differences in AX and AY are 
probably related to the extraction technique as discussed 
in Sect. |2.1| The values reported for 9 agree with those in 
this paper. This confirms the need to extrapolate for the 
in-transit orbits, in particular for 8, which earlier we con- 
cluded was the most important decorrelation parameter for 
HD 189733. 



5 SUMMARY AND DISCUSSION 

We report a re-analysis of transmission spectra observed 
with NICMOS for HD 189733, GJ-436 and XO-1. We used a 
linear decorrelation model, similar to that used in other HST 



© 2002 RAS, MNRAS 000,[T]{T7] 



A new look at NICMOS transmission spectroscopy 15 



1 I 



1.5 

/mi 



Figure 27. Transmission spectrum of XO-1 generated by deter- 
mining the planet radius from the decorrelated light curves p for 
the 18 wavelengths channels. Out results are the red points, and 
those from T10 are shown in grey for comparison. The variation 
is not caused by the atmosphere of the planet, but rather the sys- 
tematics in the data and the decorrelation method used to correct 
for these. 



analyses (e.g. |Pont etliLpO()8| S08, P09, T10), where p 



sible, to decorrelate the light curves as they all suffer from 
considerable systematics. We used a residual permutation 
algorithm in an effort to obtain more realistic uncertainties 
than in the literature, although this does not account for 
uncertainties which arise from orbit-to-orbit flux offsets. 

For HD 189733, we find a similar transmission spectrum 
to S08 when following their procedure as closely as possi- 
ble. However, using residual permutation, the uncertainties 
we calculate are often significantly larger, thus reducing the 
significance of the spectral features. We show that signifi- 
cant systematic noise remains in the in-transit orbits, which 
is not fully removed by the linear decorrelation model, or 
the channel-to-channel correction. The shape and amplitude 
of the spectral features are also strongly dependent on the 
choice of orbits, decorrelation parameters, and the nature of 
the model assumed for the baseline function. Given there is 
no physical reason to assume the baseline flux should follow 
a linear function of a particular set of optical state param- 
eters, there is no reason to prefer one model over another. 
We take this as evidence that the linear decorrelation model 
used is not a robust method to remove systematic effects 
from the light curves, and that the baseline function in the 
case of NICMOS grism data cannot reliably be given by a 
linear function of the optical state vectors from S08. 

Fig. [28] shows a plot of the overall transmission spec- 
trum of HD 189733, including the ACS optical data from 



Pont et al. (20081, the NICMOS photometric measurements 



from Sing et al. (20091, and the reduction using only orbits 



2 and 4 from this paper. This shows that we can interpret 
the NICMOS transmission spectrum as consistent with the 
optical haze, if we adopt this reduction. We emphasise that 
our analysis does not rule out the presence of molecules in 
the NIR transmission spectrum of HD 189733, but merely 
shows that the detection of molecular species is highly de- 
pendent on the decorrelation method used. 

For GJ-436 we find an amplitude of about 0.07 % ab- 



sorption in the transmission spectrum, which is similar to 
that seen in P09; however, the overall levels are inconsis- 
tent. We also find the transmission spectrum is inconsistent 
with the transit depth of the white light curve. The latter is 
much more robust given that the integrated white light curve 
suffers from fewer systemics. The transmission spectrum is 
consistent with the white light curve when only decorrelat- 
ing using the orbital phase parameters. This suggests that 
decorrelating with more parameters, particularly those with 
offsets in the parameter values, induces spurious offsets in 
the transmission spectrum. Generally, we conclude that it is 
unstable to decorrelate HST data using decorrelation param- 
eters that contain offsets between each orbit, as these offsets 
in easily be transferred to spurious offsets in flux levels of 
the in-transit orbits. The fact that we can produce a vari- 
able but unphysical transmission spectrum with the linear 
decorrelation model, suggests that we cannot extract use- 
ful spectral information below about 0.1% using NICMOS 
transmission spectroscopy, at least using current methods. 

For XO-1 we are unable to reproduce the results from 
T10, because the decorrelation parameters we measure re- 
quire an extrapolation for the in-transit orbits to determine 
the baseline function, and introduce large offsets in the flux 
levels for different orbits. The best transmission spectrum 
that we can produce show uncertainties much larger than 
those reported in T10. We do not have enough informa- 
tion to identify the source of the discrepancy, but we ex- 
pect that their uncertainties are underestimated, in a similar 
way to HD 189733. This is all the more puzzling given that 



Burke et al. (2010) also produce decorrelation parameters 



that would require extrapolation for the in-transit orbits. 
We note the importance of reporting in detail the methods 
used to correct systematics when they dominate the error 
budget. 

Our re-analysis of all three datasets shows that the sys- 
tematics in NICMOS transmission spectra cannot reliably 
be corrected for - at the level needed to detect molecular ab- 
sorption in hot Jupiters - using the multi-linear decorrelation 
techniques described in S08 and repeatedly used since. In- 
deed, there is no physical reason why the systematics should 
be described by a linear combination of the selected state 
parameters. In the absence of a better systematics removal 
technique, there appears to be a ~0.1 % floor below which 
real absorption variation cannot be distinguished from sys- 
tematics with NICMOS. All claimed detections lie near this 
level, regardless of object (Fig. [I]). Whilst this could be coin- 
cidental, more evidence is required to support these results, 
particularly when |Pont et al.] ( |2008[ ) and |Sing et al.| ( |2009 1 
suggest otherwise for HD 189733. 
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Figure 28. Overall transmission spectrum of HD 189733, including the ACS optical data from [Pont et ah] | |2008| green circles), the 
NICMOS photometric measurements from |Sing et al.| l |2009| blue squares), and a re-reduction of G206 NICMOS grism data using only 
orbits 2 and 4 from this paper (red circles) . The dotted-dashed line is the haze given in |Lecavelier Pes Etangs et al.| ([2008 ) , which the 
|Pont et al!] ( |2008| l and |Sing et al.| j2009[ | results are consistent with. The NICMOS grism data is consistent with the haze using this 
reduction. We emphasise that this case is in no way special, rather we want to illustrate that the NICMOS transmission spectroscopy 
data may still be consistent with the haze interpretation rather than the reported molecular species, given the difficulty in propagating 
the effects of systemics to the final uncertainties. 
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